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ABSTRACT 


Temporal  wind  variability  is  examined  from  the  standpoint  of  using 
wind  speed  as  the  independent  variable  since  correlations  of  wind  varia¬ 
bilities  with  height  show  no  significant  values. 

Data  from  ballistic  tests  made  at  Ft  Huachuca,  Arizona,  are  analyzed 
to  formulate  equations  for  variabilities  of  zone  and  ballistic  winds 
for  time  intervals  of  2  to  8  hours.  The  equations  involve  not  only  the 
square  root  of  the  time  (similar  to  earlier  investigations),  but  also 
the  mean  wind  speed. 

Analysis  of  the  Meppen,  West  Germany,  data  is  added  to  extend  the 
time  lag  to  1  hour  and  to  confirm  the  validity  of  the  formula  for  Western 
Europe. 
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VARIABILITY  OF  BALLISTIC  WINDS 


INTRODUCTION 

Over  the  past  20  years,  more  than  100  papers  have  been  published 
on  temporal  wind  variability.  Time  periods  have  ranged  from  seconds  to 
days,  and  levels  from  the  lowest  boundary  layer  to  the  mesosphere  have  been 
investigated.  A  strong  motivation  is  required  to  add  to  this  already 
formidable  documentation  so  that  previous  efforts  are  not  duplicated  nor 
earlier  works  merely  confirmed. 

If  one  examines  the  published  reports,  it  is  found  that  they  general¬ 
ly  follow  a  given  procedure.  Wind  differences  at  various  time  lags  are 
calculated  for  a  given  level ,  either  at  a  constant  pressure  surface 
(700  millibars  (mb) ,  500  mb,  etc.)  or  a  constant  height  (1000  meters, 

3000  meters,  etc.). 

The  dependence  of  temporal  wind  variability  upon  height  is  prevalent 
in  almost  all  meteorological  literature,  starting  with  Durst's  work  in 
1954  1  where  o,  the  root-mean-square  variation  of  the  wind  with  time, 
was  given  as  a  function  of  altitude.  Others  such  as  Arnold  and  Bellucci,2 
and  Bellucci  and  Gabriel, 3  investigating  the  lowest  30,000  feet,  showed 
a  constant  variability  dependent  only  on  the  square  root  at  the  time. 

However,  most  investigators  found  that  the  variability  increased  from  the 
lowest  levels  to  about  30,000  feet,  then  decreased  upwards.  Due  to  the 
scarcity  of  data,  variability  estimates  are  less  reliable  above  80,000  feet 
especially  for  time  lags  of  less  than  2  hours.  A  typical  chart  is  shown  in 
Fig.  1,  representing  a  composite  of  data  from  several  investigations 
(see  references  4, 5,6, 7, 8  &  9).  The  variability  maximum  at  30,000 
feet  is  clearly  shown  for  all  time  lags. 

While  all  investigations  point  to  a  maximum  of  wind  variability 
associated  with  the  level  of  the  maximum  wind,  there  are  no  reports 
describing  the  relationship  of  wind  variability  to  wind  speed. 

Mr.  A.  G.  Matthewman,  in  a  Meteorological  Working  Paper  (MWP/2) 

for  NATO  Army  Armaments  Group  AC/252  (Panel  XII),  suggests  a  slightly 

different  approach.  His  proposed  variability  formula  for  Western  Europe  is 


at  =  0.2a  t0*5  1/2  hr  <  t  <  24  hrs  (1) 


where  a  is  the  climatological  standard  deviation  (in  knots)  at  a 
particular  level. 


The  values  of  o  are  interpolated  from  U.  S.  Navy  charts10 
for  given  pressure  levels,  i.e.,  850  mb,  700  mb,  etc.  Here  a  multi¬ 
plicative  factor  is  introduced  to  account  for  the  nonconstancy  of  the 
variability  with  height.  In  essence,  the  factors  are  associated  with 
fixed  pressure  levels  above  the  surface, still  using  height  as  the  in¬ 
dependent  variable.  While  the  formula  is  an  improvement  over  the  earlier 
formula  of  Arnold  and  Bellucci  (loc.cit.). 


a  =  3.5  t0-5  1<  t  <12  hrs,  (2) 


the  physical  reasoning  behind  the  use  of  altitude  is  absent. 


DISCUSSION 

It  seems  reasonable  to  investigate  the  variability  problems  by 
considering  all  the  pertinent  parameters,  i.e.,  wind  variation,  wind 
speed,  height,  and  time  lag  to  determine  for  any  time  interval  the  para¬ 
meter  that  correlates  best  with  wind  variation.  Thus,  an  equation  of 
the  form 


AV  =  A  +  Bv  +  Ch  (3) 


where  AV  is  the  wind  difference  between  any  two  observations,  V  is  the 
mean  wind,  h  is  the  height  (A,8,C,  are  constant)  was  considered  a  more 
accurate  equation  for  the  time  variability  of  winds. 

For  this  study,  the  vector  wind  difference  between  two  observations 

v  and  y  was  used  as  a  measure  of  wind  variability, 

i  ?. 

6V  =  Vo  -  V  (4) 

2  1 

For  any  given  time-lag, the  mean  absolute  values  of  all  such  pairs 
are  computed  tc  give  the  variability  of  the  wind,  w.  Expressed 
mathematically. 


W  =  J£|  a  V|  (5) 

where  n  is  the  number  of  observations. 
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Mercator  Grid  (Zone  : 
grid  numbers  omitted 


While  w,  the  statistic  used,  is  somewhat  smaller  than  the  root- 
mean-square  value  of  the  wind  variations  used  by  many  investigators, 
tests  run  on  several  samples  of  data  show  insignificant  differences 
in  the  two  results. 

A  question  now  arises:  What  wind  speed  should  be  associated  with 
a  given  av,  since  a  pair  of  observations  taken  several  hours  apart  is 
involved.  Duplicate  calculations  were  made  using  1)  AVi,  the  earlier 
wind  speed,  and  2)  the  arithmetic  mean  of  the  two  wind  speeds,  l/2(Vi  +  V2  ). 
Variability  results  were  almost  identical  for  both  sets  of  computations. 

Data  for  the  investigation  were  taken  from  the  Artillery  Meteor¬ 
ological  Study, 11  a  field  study  made  in  the  Fort  Huachuca,  Arizona,  area 
during  January  and  February  1965.  A  meso  network  of  12  rawin  stations 
was  established  to  obtain  synoptic  soundings  in  conjunction  with  the 
firing  of  two  8-inch  howitzers.  Simultaneous  balloon  flights  were 
made  at  2-hour  intervals  from  0600  LST  to  1400  LST  on  16  test  days, 
a  total  of  960  soundings.  Data  were  reduced  by  computer  after  careful 
editing  to  remove  reading  and  punching  errors  and  gross  errors  in  the 
angular  recorder  print-out. 12  The  stations  and  their  locations  are 
shown  in  Fig.  2. 

Since  the  field  tests  were  intended  for  artillery  purposes,  winds 
were  computed  for  the  artillery  zones,  as  shown  below  in  Fig.  3.  The 
mean  wind  for  the  layer  was  determined  and  that  yalue  attributed  to  the 
mid-point  of  the  zone,  i.e.,  for  Zone  1  the  wind  is  assigned  to  350  meters. 
Similarly,  computations  were  performed  using  standard  artillery  weighting 
factors **  in  order  to  determine  the  variability  of  ballistic  winds. 
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Figure  3.  Artillery  Zone  Structure. 
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Results 

Returning  to  the  general  variability  equation  given  earlier,  Eq.  (3), 
AV  *  A  +  BV  +  Ch,  coefficients  A,  B,  and  L  are  determined  by  standard 
linear  regression  techniques, 17  and  their  relative  size  indicates  the 
degree  of  dependence  of  AV  on  V  or  H.  Table  1  presents  the  values  of  the 
coefficients  for  the  indicated  time  lags,  along  with  the  number  of  pairs 
used  in  the  computations. 


TABLE  1 

Regression  Coefficients:  Unweighted  Zone  Winds 
AV  =  A  +  BV  +  CH 


t 

(hours) 

n 

A 

(knots) 

B 

(knots) 

C 

(knots) 

2 

8586 

3.75 

0.12 

-0.00015 

4 

6303 

5.49 

0.16 

-0.00022 

6 

4160 

6.74 

0.22 

-0.00030 

8 

2019 

7.67 

0.27 

-0.00043 

It  is  obvious  that  the  coefficients  vary  with  time  lag,  and 
furthermore,  the  variation  is  not  linear.  To  determine  the  functional 
relationship,  plots  of  av  against  the  coefficients  were  constructed 
on  log-log  paper.  The  resultant  slope  indicated  that  the  exponent  of  "t" 
should  be  0.5.  The  equations  for  each  time  lag  are  shown  below: 


Variability 

of  Zone 

!  Winds 

At  = 

2 

hrs 

Av  = 

(2.65 

+ 

0.08V  - 

■  0.0001 

H)t°< 

,5 

At  = 

4 

hrs 

Ay  = 

(2.75 

+ 

0.08V  - 

•  0.0001 

H)t°- 

,5 

At  = 

6 

hrs 

Av  - 

(2.75 

+ 

0.09V  - 

•  0.0001 

H)t°< 

,5 

At  = 

8 

hrs 

Av  = 

(2.72 

+ 

0.09V  - 

■  0.0002 

H)t°  - 

.5 

6 


Combining  these  equations  into  a  single  equation  we  obtain: 


AV  =(2.7  +  0.08V  -  0.  0001  H)t°*5  2  <  t  <  8  hr  (6) 


Equation  (6)  makes  it  very  clear  that  there  is  virtually  no  depend¬ 
ence  of  wind  variation  with  height.  This  is  emphasized  when  the  correla¬ 
tion  coefficients  are  examined  as  indicated  in  Table  II. 


TABLE  II 


Partial  Correlation  Coefficients  (Zone  Winds) 


— 

t 

(hours) 

*12 

R13 

R2  3 

2 

0.42 

0.08 

0.44 

4 

0.46 

0.09 

0.45 

6 

0.50 

0.09 

0.45 

8 

0.56 

0.08 

0.46 

R12  =  correlation  coefficient  between  av  and  H 

R13  =  partial  correlation  coefficient  between  av  and  H 

R23  =  partial  correlation  coefficient  between  av  and  H 


Note  that  the  wind  variability  correlates  well  with  the  wind  speed 
(R  -  0.5),  but  not  at  all  with  height  (R<0.1),  while  the  wind  speed 
itself  correlates  with  the  height.  It  is  this  last  correlation  that 
causes  an  apparent  relationship  of  wind  variability  with  height,  which 
is  noted  so  often  in  the  literature. 

Ballistic  Winds 


When  the  zone  winds  are  weighted  with  the  standard  weighting  factors  12 
to  give  an  integrated  wind,  a  smoothing  results  which  should  tend  to  reduce 
the  wind  variability.  This  is  indeed  the  case  as  noted  below. 


Van  ability  of  Ballistic  Winds 


At  =  2  hrs 

4VB 

=  (2.49  +  0.07Vb-  0.00016  H)t0-5 

At  =  4  hrs 

avb 

=  (2.35  +  0.08Vb-  0.00014  H)t0*5 

At  =  6  hrs 

avb 

=  (2.38  +  0.09Vb-  0.00013  H)t0-5 

At  =  8  hrs 

avb 

=  (2.36  +  0.09Vb-  0.00013  H)t0-5 

where  aVjj  is  the  variability  of  ballistic  winds,  and  Vg  is  the 
ballistic  wind  speed. 

Again,  combining  to  form  a  single  equation: 


AVb  =  (2.4  +  0.08Vb  -  0.0001  H)t0-5  2  <  t  <  8  hrs  (7) 

Note  that  V  from  Eq.(6)  and  Vg  from  Eq - ( 7)  are  different.  In  general, 
the  ballistic  wind  is  about  80%  of  the  unweighted  wind  for  the  comparable 
zone.  In  this  study  the  mean  values  of  the  zone  and  ballistic  winds  were 
26.5  knots  and  21  knots  respectively,  conforming  well  to  the  ratio 
established  from  other  studies. 

The  correlation  coefficients  are  not  significantly  different  from 
those  of  Table  II,  as  noted  from  Table  III  below: 


TABLE  III 

First-Order  Correlation  Coefficients  (Ballistic  Winds) 


t 

(hours) 

Ri?. 

R13 

,j»  \ 

to 

1 

1 

i 

1 

i 

2 

0.34 

0.08 

0.49 

4 

0.46 

0.02 

0.49 

6 

0.53 

0.06 

0.49 

8 

- 

0.57 

0.08 

0.49 

The  almost  total  lack  of  correlation  between  wind  variability  and 
height  would  lend  credence  to  the  dropping  of  the  "H-term"  in 
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the  variability  equations.  In  the  case  of  application  to  artillery 
where  all  maximum  ordinates  are  less  than  20  km  and  most  below  10  km, 
there  is  additional  justification.  The  simplified  variability 
equations  become; 


AV  =  (2.7  +  0.08V)t°  *5  2  <  t  <  8  hrs  (8a) 

aV  =  (2.4  +  0.08V..  )t0*5  2  <  t  <  8  hrs  (8b) 

B  o 

It  will  be  noted  that  the  older  Eq.  (2)  on  wind  variability, 
which  was  an  average  value  of  many  observations,  may  be  considered  as 
the  wind  variability  associated  with  a  10-knot  wind,  which  is  a  good 
estimate  of  the  mean  wind  observed  during  the  investigation. 

Since  meteorological  data  from  the  Fort  Huachuca  tests  covered  a 
period  of  2  to  8  hours,  extension  of  the  period  of  validity  of  Eqs. 

(8a)  and  (8b)  from  such  data  is  not  possible.  Some  references,  however, 
can  be  made  from  general  considerations.  It  is  obvious  that  Eqs.  (8a) 
and  (8b)  cannot  be  used  for  long  time  intervals,  since  the  wind  vari¬ 
ability  would  take  on  impossibly  large  values.  A  different  form  of 
variability  equation  must  be  employed  that  follows  from  the  standard 
vector  deviation  of  the  difference  of  two  vectors: 


al-2  =(°12  +  °22  +  °l02rl,2)0,5  (9) 


where  oj,  a?  are  standard  deviations  of  two  vectors,  and  r  is 
the  correlation  coefficient.  1  z 

In  the  case  of  wind  variability  with  time: 


°l-2  ~  0 


Oj  =  02  =  O 


For  very  long  time  periods,  r  0,  hence: 


=  /  2  a 
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Thus  the  variability  approaches  *^2  times  the  standard  deviation  of 
the  wind  as  a  limit. 

The  Meppen  trials14  use  radiosonde  data  with  longer  time  separations 
between  flights  -  up  to  24  hours.  An  analysis  of  those  data  are  contained 
in  the  Addendum  to  this  report. 

At  the  other  end,  investigations  of  R.  Bellucci15  have  shown  th„t 
tor  very  short  time  intervals  (less  than  one  hour)  wind  variabilities  do 
not  follow  the  "sguare  root"  law.  Thus,  Eqs.  (8)  should  not  be  used  for 
t  <1.  Since  Eq.  (2)  of  Arnold  and  Bellucci  indicates  a  lower  limit  of 
1  hour,  the  lower  limit  of  our  Eqs.  (8)  is  tentatively  set  at  1  hour 
pending  verification  from  Meppen  Analysis. 

The  equations  now  become: 


AV  =  (2.7  +  0.08V)t°*5  1  <  t  <  8  hrs  (12a) 

AV  =  (2.4  +  0.08V  )t°* 5  1  <  t  <  8  hrs  (12b) 

D  B 

CONCLUSIONS 

In  this  report,  a  new  equation  for  wind  variability  has  been  derived 
for  periods  between  1  hour  and  8  hours,  which  depends  not  only  on  the  time 
interval  between  observations,  but  also  on  the  wind  speed  itself.  It  would 
be  desirable  to  have  a  complete  formulation  for  all  periods,  but  as  noted, 
the  time  lag  becomes  less  important  at  short  as  well  as  at  very  long  periods. 


Equation  (10)  should  be  made  applicable  for  all  periods,  but  the 
difficulty  in  determining  the  correlation  coefficient  makes  it  somewhat 
unwieldly.  An  assumption  is  generally  made  that  the  correlation  coeffi¬ 
cient  is  exponential,  viz.: 


r  =  e-ctt  (13) 

This  is  probably  an  oversimplification  since  investigations  1,9,16 
show  that  r  varies  with  geography,  season,  altitude,  time  lag,  etc.  In 
addition,  for  time  lags  larger  than  36  hours,  there  is  a  systematic  de¬ 
parture  of  "r"  from  Eq .(13).  A  re-examination  might  be  indicated  using 
wind  speed  as  the  independent  variable. 

Variability  studies  for  short  time  lags  are  difficult  because  of 
the  requirement  for  multiple  ground  equipments  and  the  inherent  measure¬ 
ment  errors  that  effectively  mask  the  small  variability  associated  with 
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periods  of  a  few  minutes.  A  study  of  wind  variability  in  the  2-  to  10- 
minute  range  conducted  by  R.  Bellucci  (op.cit.)  using  pibal  balloons, 
showed  that  wind  variabilities  were  approximately  25%  of  the  wind  speed 
and  essentially  independent  of  time  in  the  2* to  10-minute  range. 

The  relative  unimportance  of  time  in  short-period  wind  variability 
may  be  noted  from  examination  of  the  work  of  Plagge  and  Smith  (op.cit.). 

The  values  of  the  time  variability  of  wind  for  one-quarter  hour  and  one- 
half  hour  are  not  significantly  different,  and  definitely  do  not  vary  as  the 
square  root  of  time  lags. 

Thus,  a  complete  explicit  formulation  of  wind  variability  in  terms 
of  wind  speed  and  time  lag  would  be  complex,  involving  at  least  three 
terms.  The  first  term  would  apply  to  short  time  periods,  another  to 
very  long  periods,  and  the  third  term  would  be  similar  to  our  Eqs.  (8). 
Beyond  its  period  of  validity,  the  value  of  each  term  would  be  almost 
negligible  compared  to  the  applicable  terms,  hence  its  contribution  is 
felt  only  when  required. 

Symbolically 

AW  =  S  +  M  +  L  (14) 

where  the  value  of  S  0  when  t  >  1  hour, 

and  the  value  of  M  •>  0  when  t  <  1  hour  or  t  >  8  hours, 

and  the  value  of  L  -*■  0  when  t  <  8  hours. 

The  limits  of  "t”  are  meant  to  be  illustrative  rather  vhan  definitive. 

Considerably  more  data,  especially  at  the  shorter  time  intervals,  are 
required  before  any  reliable  relationships  can  be  established. 


ADDENDUM 


Analysis  of  Meppen  Data 


The  Meppen  data  from  the  Summerwind  Program  were  analyzed  following 
the  procedures  outlined  earlier  in  the  report,  to  the  extent  that  the 
data  would  allow.  In  that  program,  nine  upper  air  stations  were  estab¬ 
lished  with  soundings  made  at  intervals  as  short  as  1  hour.  While  the 
instrumentation  for  measuring  the  winds  was,  in  most  cases,  different 
from  that  employed  at  Fort  Huachuca,  the  care  with  which  the  data  were 
taken  makes  the  inclusion  of  this  work  worthwhile. 

Results 

The  regression  coefficients  for  ballistic  winds,  with  the  number 
of  paired  values,  are  shown  in  Table  IV.  It  will  be  noted  that  the 
number  of  cases  is  far  less  than  in  the  Fort  Huachuca  Study,  and  that 
there  is  a  sharp  decrease  between  6  and  8  hours.  The  number  of  data  pairs 
remains  small  until  the  24-hour  interval,  with  the  attendant  loss  in 
significance  of  the  results  between  8  and  15  hours. 

TABLE  IV 

Ballistic  Winds  (Upper) 

Regression  Coeffi cents 
AV  =  A  +  BV  +  CH 


t 

!  n 

A 

B 

C 

(hours) 

f 

i 

(knots) 

1 

'  2135 

2.45 

0.085 

1 

1  2 

1422 

3.11 

0.081 

0.00028 

!  3 

1  2077 

3.16 

0.132 

0.00021 

4 

.  1925 

4.13 

0.119 

0.00018 

6 

.  1314 

3.75 

0.200 

0.00040 

8 

160 

3.21 

0.159 

0.00031 

10 

1  417 

3.05 

0.294 

0.00001 

12 

|  255 

1.22 

0.455 

-0.00017 

15 

228 

0.68 

0.459 

0.00010 

24 

1357 

-1.14 

0.471 

-0.00015 

Proceeding  as  before,  we  form  equations  for  each  time  lag  by 
dividing  out  t° -f  The  results  are  shown  on  the  following  page: 
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Variability  of  Ballistic  Winds 


t  =  1  hr 

AVB 

=  (2.45  +  0.09  VB  +  0.00013  H)t° -5 

t  =  2  hrs 

aVb 

=  (2.20  +  0.06  VB  +  0.00019  H)t0*5 

t  =  3  hrs 

aVb 

=  (1.83  +  0.08  VB  +  0.00012  H)t0*5 

t  s  4  hrs 

aVb 

=  (2.06  +  0.06  VB  +  0.00010  H)t0*5 

t  =  6  hrs 

aVb 

=  (1.53  +  0.08  VB  +  0.00002  H)t0*5 

t  =  8  hrs 

AVB 

=  (1.14  +  0.06  VB  +  0.00011  H)t0,5 

t  =  10  hrs 

avb 

*  (0.97  +  0.09  VB  +  0.00003  H)t°-S 

t  =  12  hrs 

aVb 

=  (0.35  +  0.13  VB  -  0.00005  H)t0-5 

t  =  15  hrs 

AVb 

=  (0.17  +  0.12  VB  +  0.00025  H)t0*5 

t  =  24  hrs 

AVb 

=  (-0.2  +  0.10  VB  -  0.00003  H)t0*5 

Comparison  of  the  above  results  and  those  given  earlier  shows  ex¬ 
cellent  agreement  up  to  6  hours  and  only  fair  agreeement  at  8  hours.  For 
long  periods,  Eq.(7)  is  obviously  invalid.  The  discrepancy  at  8  hours 
is  strange  since  a  similar  time  lag  at  Fort  Huachuca  gives  results  as 
predicted.  However,  there  is  one  order  of  magnitude  difference  in  the 
number  of  cases.  It  may  well  be  that  160  cases  are  not  a  representative 
sample  for  the  statistic. 

Weighting  the  various  time  lags  to  form  one  equation,  we  obtain: 


AVg  =  (2.0  +  0.07  Vg  +  0.0001  H)t0*5  1  <  t  <  8  hr  (15) 

Finally,  simplifying  by  ignoring  the  height  term, 

AVg  =  (2.0  +  0.07  Vg)t0-5  1  <  t  <  8  hrs  (16) 

As  postulated  earlier,  we  have  confirmed  that  the  Eqs.  (8)  are  valid 
for  i  hour.  The  upper  limit,  8  hours,  is  still  within  acceptable  con¬ 
fidence  limits,  the  standard  error  of  estimate  being  4.7  knots. 
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The  concept  that  wind  variability  for  longer  periods  does  not  follow 
the  "t0*5"  law  is  very  evident  from  the  Meppen  data  since  there  is  scarcely 
any  difference  in  the  12-  and  24-hour  wind  variabilities.  If  additional 
data  were  available,  the  variabilities  would  continue  to  increase  very 
slowly  since  the  limiting  value,  aJl  ,  has  not  yet  been  reached. 

The  value  of  the  constant  “V  ~|0.2|  for  the  two  longest  time 
lags  (15-  and  24-hours)  is  insignificant  compared  to  the  "BV"  term. 

It  is  noteworthy  also  that  the  coefficient  of  V  is  relatively  invariant, 
approximately  one- tenth  of  the  mean  wind  speed.  This  contrasts  strongly 
with  the  variability  for  very  short  periods  (~  few  minutes),  where  the 
dependence  is  of  the  order  of  one-fourth  or  one-fifth  the  wind  speed. 

SUMMARY 


The  Meppen  data  confirms  the  validity  of  the  proposed  wind  vari¬ 
ability  formula  for  the  period  of  2  to  8  hours,  and  extends  the  lower 
limit  to  1  hour.  The  data  likewise  show  that  the  formula  cannot  be 
used  for  time  lags  greater  than  12  hours. 
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